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O I 1 Abstract 

Oh' We examine the dark matter phenomenology in supersymmetric light higgs boson scenarios, 

adapting nonuniversal Higgs masses at the gauge coupling unification scale. The correct relic 
density is obtained mostly through the annihilation into a pseudoscalar A, which gives high 
values for the self-annihilation cross-section at present times. Our analysis shows that most 
part of the A pole region can produce detectable gamma-rays and antiproton signals, and still 
J^\ . be compatible with with recent direct detection data from XENON100 and CDMS-II. 
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1 Introduction 

Low energy super symmetry (SUSY) pQ, has emerged as the most promising candidate for physics 
beyond the Standard Model (SM). The minimal super symmetric standard model (MSSM) has been 
attractive for its explanations towards several phenomena such as (i) the Higgs mass hierarchy 
problem in the SM, (ii) electroweak (EW) symmetry breaking via radiative corrections and (iii) the 
unification of the SM gauge couplings at the grand unification scale (Mqut — 2 x 10 16 GeV). 
Assuming R-parity conservation, a specially attractive feature of the MSSM is the presence of the 
stable, weakly interacting lightest super symmetric particle (LSP) in the form of lightest neutralino 
(x?) [2] which turns out to be a good candidate for the observed cold dark matter (CDM) in the 
universe. 

However, the MSSM has a large number of soft breaking parameters, which make it less predic- 
tive in its most general forms. Most of these parameters can be elliminated by considering a definite 
SUSY breaking framework at the high scale. The simplest model in this regard, which assumes 
gravity as the mediator for SUSY breaking - the minimal supergravity (mSUGRA) [3] model - has 
only five free parameters. These are the common gaugino and scalar mass parameters m\/2 and 
mo, the common tri-linear coupling parameter Aq, all given at the gauge coupling unification scale, 
the ratio of Higgs vacuum expectation values at the electroweak scale, namely tan /3, and the sign 
of the Higgsino mixing parameter /i. The lightest neutralino can be expressed in terms of its gauge 
eigenstates, the Bino, the Wino and the two Higgsinos as 

X ° 1 = N 11 B + N 12 W + N 13 H D + N 14 H U . (1) 

The model is primarily associated with a Bino-dominated LSP that leads to over-abundance of 
dark matter. Reduction of the relic density to satisfy the WMAP data can be achieved if i) there is 
coannihilation of the LSP with another sfermion (usually stau f\ or rarely t% ) that has a mass close 
to the mass of the LSP |3], ii) there is appropriate mixture of Bino and Higgsinos in the composition 
of the LSP so that there may be coannihalating charginos in the LSP-Xi annihilation [5l [6] (the 
XiXi^ coupling is also enhanced), the so called focus point (FP) [7j /hyperbolic branch (HB) [8] 
region (this kind of coannihilation becomes significant when the LSP is dominantly Higgsino like), 
iii) the LSP is sufficiently small in mass and sfermions are light so that light sfermion exchange may 
enhance the LSP-LSP annihilation rates p, or iv) there is a possibility of having s-channel Higgs 
exchanges occuring via the CP-odd Higgs boson A or via CP-even heavy (light) Higgs bosons H (h) 
leading to the "funnel region" of the dark matter satisfying zone in the m\/2 — mo plane |10|. II 1 j . 
This simple model is highly predictive because of its small number of unknown parameters. 
However, a very severe restriction appears on the {rfi\ji — mo) plane of mSUGRA from the LEP2 
bound on lightest Higgs boson mass (roji) [12]; which reads 

m h > 114.4 GeV. (2) 



lr This region known as bulk annihilation region, though highly constrained by the LEP2 limits on the Higgs boson 
mass, could be revived with non- vanishing Ao parameter [S]. 



In order to satisfy the limit m^ >114 GeV within the framework of mSUGRA, one requires rea- 
sonably large values for the top squark mass (S> Mz) which implies some fine-tuning of SUSY 
parameters to obtain the correct value of Mz- However, the limit of Eq.Q is related to the 
principle Higgs search channel in the LEP2 that involves production of the Higgs boson via the 
e + e~ — > Zh process, followed by decays of the Higgs boson. The production process principally 
depends on the coupling gzzh- Now, in the MSSM, the coupling gzzh is not identical with the 
respective coupling in the SM, it can rather be written as sin(/3 — a)g^ h (For reviews see |13||l4j). 
Here, gzz h represents the SM Higgs coupling to the Z boson. The prefactor sin(/3 — a) appears 
due to the mixing of the CP-even Higgs bosons, where a is the mixing angle. In supersymmetric 
theories the limit of Eq.([2]) refers to the scenarios where sin(/3 — a) ~ 1. This corresponds to the 
so-called decoupling zone for the Higgs boson where one finds ra\ » M§. In fact this is true for 
most of the parameter points in mSUGRA. However, in the general MSSM the above prefactor 
may become much smaller. This opens up a supersymmetric parameter space associated with a 
Higgs mass smaller than the limit of Eq.([2]) (For details see |14|). This is the non decoupling zone 
for the Higgs boson, where all Higgs masses may have comparable values mh ~ mu ~ tua ~ Mz 
|15 1 116 1 [T7 1 118 1 ITU 1 120 1 f2T | 122] . An important virtue is that with this new window, a light Higgs mass 
with rah ~ (95 — 100) GeV could explain the slight excess of events of 2.3a statistical significance 
|15j as reported by the LEP2 experiments [12]. Within this non decoupling zone it is possible to 
have a Higgs boson lighter than the Z boson [161 117] . In this regime, the Higgs mass does not 
require large radiative corrections, thus could evade the little hierarchy problem pL8] . 

This particular region, which has recently been named as ^ Light Higgs Boson 1 scenario or LHS 
[171 H2]> is a iso interesting as it provides definite predictions. Since all Higgs states are light 
enough, they can be produced at the Large Hadron Collider (LHC) [T7J. Interestingly, this region 
can also accommodate the WMAP data for the dark matter relic density [19] , However, realization 
of the light Higgs boson scenario is not possible in the most constrained model like mSUGRA. 
Scenarios like supergravity - inspired nonuniversal Higgs mass models (NUHM), where the Higgs 
mass square parameters may have different values at the scale Mqut (w?h (Mgut),™<h {^gut)) 
from the other scalars |231 124]. may be useful in this regard. Similarly, this particular SUSY zone 
can be realized in a particular nonuniversal scalar mass model where third generation squarks and 
the two Higgs doublets may have vanishing values at the GUT scale |25j . 

In this work, adapting a NUHM framework, we try to explore this particular parameter space in 
the light of recent dark matter detection experiments. Moreover, our focus will be to compute the 
indirect detection rate, particularly where the mass of the lightest neutralino and Higgs bosons are 
reasonably small. We will try to find out the parameter space where the WMAP data is satisfied 
via the resonant annihilation into a CP-odd Higgs boson. This is interesting as then neutralino pair 
annihilations in our galaxy could produce detectable signals both in the antiproton and the gamma- 
ray channels. In fact, the part of the MSSM parameter space with light neutralino masses which is 
otherwise forbidden in models like mSUGRA, can completely be explored in the near future. The 
other WMAP - satisfying zones, as mentioned, cannot produce such large anti-proton/gamma ray 
signals (except the FP/HB zone where the Bino dominated LSP with sufficient Higgsino component 



could lead to large signals). Denning a measure for these observables, we examine the detection 
prospects for the parameter space in the ongoing as well as upcoming experiments. 

Now, it is not surprising that the neutralino - nucleus spin-independent scattering cross-section 
a SI , which is dominated by i-channel Higgs exchange (h,H) could be very large in the present 
context. The cross-section reaches its maximum value in the non decoupling zone because of 
smaller mass values of the CP-even Higgs bosons. In fact this particular region may fall within the 
range of the exclusion limits as given by CDMS-II [26J or XENON100 [27J. However there exist 
considerable uncertainties in the computation of the above cross-section as well as the relevant 
experimental limits. The principle sources are related to (i) nuclear form factors, (ii) the WIMP 
- nucleon scattering cross-section (mainly concerning the strange quark content of the nucleon), 
(iii) the local dark matter density and (iv) the escape velocity of the WIMP. Considering all these 
uncertainties and particularly exploiting the uncertainty related to the strange quark content of the 
nucleon we will see that the predicted neutralino-nucleus scattering cross-section <j SI may change 
significantly. These uncertainties render the question of exclusion of interesting regions of the 
parameter space related to dark matter phenomenology subtle to answer. 

The paper is organized as follows. In Secj2] we will primarily discuss the model and the low 
energy constraints which are relevant here. In Secj3]we introduce the spin-independent neutralino- 
nucleus scattering cross-section for the LSP. Similarly, SecJH is devoted to the indirect detection of 
the LSP, primarily to the (i) gamma-ray and the (ii) anti-proton signal. We describe our results in 
SeciSJ In SecEJwe provide estimates regarding the effect of the uncertainties in the computation of 
a SI . Finally we conclude in Seed 

2 The NUHM and low energy constraints 

In this section, we describe the model parameters and the relevant low energy constraints. As 
mentioned in the introduction, we place ourselves in a supergravity framework with nonuniversal 
Higgs masses at the GUT scale. Thus, the model can be specified by the following parameters: 

mi/2, A , m , m 2 Hu (M GUT ),m 2 HD (Mqut), sign(fx),tan j3. (3) 

Here, ra\j2 and A§ represent universal values for the gaugino masses and the trilinear coupling 
parameters, while mo denotes universal masses for all scalars except the Higgs bosons at the GUT 
scale. In general, Higgs scalars can belong to different multiplets of a grand unified gauge group, 
thus may assume different values than the other scalars in the theory (see |28j for motivated 
constructions leading to specific nonuniversal Higgs mass terms). In this spirit, we scan over the 
Higgs mass square parameters m 2 H (Mgijt) and m^ (Mqut) at the scale Mqut to produce a 
MSSM parameter space consistent with all low energy constraints and where Higgs bosons can be 
reasonably lighter. Using the Renormalization Group Equations (RGE) for the soft SUSY breaking 
parameters along with the conditions for radiative EWSB one can determine \x and ijia through 

, _ 1 2 rn% D -m 2 Hu tan 2 fi ^-^tan 2 /? 
H -~-M z + tan2/3 _ 1 + tan2/ 3_i W 
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and 

sin 2/3 = 2Bfi/{m 2 HD + m\ v + 2/x 2 + Ei + S 2 ) (5) 

where the Sj's represent the one- loop corrections |29t l30|. which become small at the scale where 
the Higgs potential Vmggs is minimized. We can approximate fi 2 ~ ~' in \ (f° r tan/3 > 5) and 
m\ = m 2 H +m 2 H + 2// 2 ~ m 2 H —in 2 H & t tree level. The parameters [i and the sparticle spectrum 
have been computed with SuSpect [31j . We consider the following collider constraints in addition 
to the bounds on sparticle masses [32] : 



Higgs boson mass limit: In the non decoupling region where the A boson becomes very 
light so that one has m^ ~ mn ~ m^ ~ Mz, the lower limit of m^ goes down to 93 
GeV or even lower. We consider that the parameter space with sin 2 (/3 — a) < 0.3 (or, 
sin(/3 — a) < 0.55), but with Higgs mass 93 < rrih < 114 is in agreement with the LEP2 
limit [12] . In practice, relaxing the constraint marginally we discern this zone with a value 
sin(/3 — a) < 0.6. Consequently, the coupling of the heavier Higgs boson to the Z boson 
(dzzH on cos(/3 — a)) becomes dominant and this makes the heavier Higgs boson SM - like 
(cos(/3 — q) ~ 1) with the same lower bound as in Eq.([2]). In summary, to obtain acceptable 
SUSY spectra with 93 < m^ < 114, in addition to the desired value for sin(/3 — a), one also 
requires m# > 114 GeV. 

On the other hand, in the decoupling region (sin(/3 — a) ~ 1 ), the limit of Eq.([2]) needs to 
be respected. However we note that there is an uncertainty of about 3 GeV in computing 
the mass of the light Higgs boson [33] • This theoretical uncertainty primarily originates from 
momentum-independent as well as momentum-dependent two-loop corrections and higher 
loop corrections from the top-stop sector. Consequently, a lower limit of 111 GeV is often 
accepted for the SUSY light Higgs boson mass. 

Br(b — > sj) constraint: The most significant contributions to b — > S'j originate from charged 
Higgs and chargino exchange diagrams in models like mSUGRA. The charged Higgs (H~-t 
loop) contribution has the same sign and comparable strength with respect to the W~-t loop 
contribution of the SM, which already saturates the experimental result. Hence, in scenarios 
where the charged Higgs mass can be very small, satisfying the b — > sj constraint requires a 
cancellation between the charged Higgs and the chargino exchange diagrams. We will choose 
Aq such that most of the parameter points, especially WMAP - compliant ones, can satisfy 
this constraint. The effect of this constraint in the light Higgs boson zone has been discussed 
in [18]. We have used the following 3<r level constraint from b — > s"f with the following limits 

2.77 x 10~ 4 < Br(b -+ s'j) < 4.33 x 10 -4 . (6) 

Br(B s — > /i + /_i _ ) constraint: Similarly, the flavor physics observable B s — > /i + /i~ may 
become very significant in this particular parameter space. The current experimental limit 



for Br{B s — > fi + /j, ) coming from CDF [35] can be written as (at 95%C.L.) 

Br(B s -> n + n~) < 5.8 x 10~ 8 , (7) 

which has recently been improved to < 4.3 x 10~ 8 at 95% C.L [36]. The estimate of B s — > 
H + H~ in the MSSM depends strongly on the mass of A-boson and on the value of tan/3. In 
particular, the neutral Higgs boson contribution scales as m~ A whereas there is a (tan/3) 6 
type of dependence. However, in the present analysis, we choose tan/3 (= 10) which makes 
this constraint less restrictive for most of the parameter space. 

• WMAP constraint : In computing the relic density constraint, we consider the following 
3cj limit of the WMAP data [27], 

0.091 < ttcDMh 2 < 0.128. (8) 

Here QcDhdh 2 is the dark matter relic density in units of the critical density and h = 0.71 ± 
0.026 is the Hubble constant in units of 100 Km s _1 Mpc -1 . We use the code micrOMEGAS 
[38] to compute the neutralino relic density. 

3 Direct Detection of Dark Matter 

Direct detection of the LSPs involves measurement of recoil energies of nuclei as they are scattered 
by WIMPs. The effective Lagrangian that describes Xi ~ Q elastic scattering at small velocities is 
given by (for reviews see |39j) 

£ = <iXi7 M 7 5 Xi%7^7 5 % + aqiXiXi^ • ( 9 ) 

The first term represents spin-dependent scattering while the second term refers to spin-independent 
scattering. Eq.Q assumes summing over both the quark generations q while the subscript i runs 
for up (i = 1) and down type (i = 2) quarks respectively. The neutralino-quark coupling coefficients 
a q and a' contain all SUSY model-dependent information. The spin-independent scattering cross- 
section of a neutralino with a target nucleus of proton number (atomic number) Z and neutron 
number A — Z (A being the mass number) is given by 

a SI = ^ [zfp + {A _ z)fn] 2 (1Q) 

m m N 

Where m r is the reduced mass defined by m r = j— — 1 , ., and rriN refers to the mass of the nucleus. 

x i 
The quantities f p and f n contain all the information of short-distance physics and nuclear partonic 

strengths. These are given by 

U,(n) _ V- f(p,(n)) ®q 2 ( p ,( n )) 






PA n ) q=u,d,s q c,b,t q 



where ff 1 ' defined as 



m p,(n)fTq = (P) ( n )\ m qQl\P, ( n )) = m q B q ■ ( 12 ) 



The quantities f^' can be evaluated using a few hadronic data [iQj. The gluon - related part 
namely f^G is given by 

/&<"» = 1 - E /^ g (n)) • (13) 

q=u,d,s 

The numerical values of f^l may be seen in [40|, I41j . We compute the neutralino-nucleon spin- 
independent scattering cross-section by using the code DarkSusy [32]. We should note here that 

(p (n)) 
the parameter /^ requires the information of the pion-nucleon sigma term a n ^ and the size 

of the SU(3) symmetry breaking-o"o as /^ oc (<j w n — o~o), so that the leading contribution in 

(jSl g oeg ag ^ (ct^jv — o"o) 2 - Clearly, uncertainty in the computation of tj^N can produce significant 

(p (n)) 

shift to the strange quark parameter frfg and consequently to the measured spin-independent 
cross-section. Recent lattice results hint towards much smaller values of the f^l ~ 0.02 [201143]. 
a value much smaller than previous estimates. Considering even larger uncertainty in ovjv one may 
assume a^N = o"o which leads to f^ = [33]. In DarkSusy [32] the above coefficient is chosen as 
j!\pa )) = 0.14. This could provide a significant change in the results of the a SI . In fact, in [2UI 144] . 

(p (n)) 

the variation in the spin-independent cross-section due to this reduced f^ has been estimated. 

We shall also compute the variation in the spin-independent cross-section with the reduced values 

(p M) 
of ftp* . Similarly, we shall comment on the other sources of uncertainties and consequently their 

effects on a and the corresponding experimental limits. 

The scalar cross-section depends on t-channel Higgs exchange (h, H) and s-channel squark 

exchange diagrams (a SI ~ — jl — ). Now, considering the strong bounds on the light squark masses, 

m H,h,q 

a SI is dominated by the exchange of CP-even Higgs bosons. In the present context, where both 
CP-even Higgs masses are quite small, a becomes enhanced. In fact, we shall see that in the 
mass range 93 < rrih < 114, the spin-independent cross-section becomes maximal and this renders 
the light Higgs boson zone highly constrained even for tan/3 = 10. Apart from the masses of the 
Higgs bosons, the cross-section depends strongly on the couplings XiXih(H) (C o a h i H \) and also 
on h(H)qq (C^^m). The latter contributions in the context of down - type fermions include tan/3 
as Cqqh ~ tan/3cos(/3 — a) and C q qH ~ tan/3sin(/3 — a). Clearly, both couplings may become quite 
large for larger values of tan/3. These large couplings, in addition to the small mass values for 
mh,H, make this light Higgs boson zone almost ruled out for larger values of tan/3. 

4 Indirect Detection of Dark Matter 

Indirect detection techniques are based on the detection of primary or secondary particles resulting 
from dark matter (DM) annihilation or decay. The observed flux of these particles is proportional 
to the annihilation rate of the DM species which in turn is, for the annihilating case, proportional 



to p 2 , the dark matter density profile. In the following, we summarize the basic ingredients of the 
formalism that we shall be employing in this work. 

4.1 Gamma-ray Detection 

4.1.1 The flux at intermediate latitudes 

The differential flux of gamma-rays generated from dark matter annihilations and coming from a 
direction forming an angle tp with respect to the galactic center (GC) is 

where N = 1, 1/2 if the annihilating particle is Majorana-like or Dirac-like respectively, (av) is 
the total self annihilation cross-section averaged over velocity for v — > 0, m x is the mass of any 
DM particle , dNf/dE is the differential spectrum of final state / into gamma-rays of energy E 
and p is the dark matter density profile. We note that in this work we only focus on prompt 
photons, ignoring potential contributions from synchrotron emission or inverse Compton scattering 
of charged particles with the interstellar medium. 

The functions -jJ- have been computed using the fragmentation of Standard Model particles 
into 7-rays according to the PYTHIA |45j Monte Carlo code. 

Now, of course, we don't observe annihilations along a line but rather within a conical region 
around some angle i/jq. It is convenient to define the dimensionless quantity 



JO/0 = — - / p 2 (l(i>))dl(iP), (15) 



1 1 

Ro Pq Jlos 

where Rq is the sun's distance from the galactic center and po the local DM density (actually, they 
can be any arbitrary normalization factors rendering J a dimensionless quantity). 
Then, we can define the average of this quantity within a solid angle A$7 as 

J(Afi,W = ir/ J(^)dfi, (16) 

and we can then rewrite the flux as 

The most common region of the sky that is examined in the literature as a source of 7-rays 
is the galactic center, since it is the region where iV-body simulations predict a maximization of 
the dark matter density distribution and, hence, the corresponding gamma-ray flux. However, the 
galactic center is a region which is known poorly : There are large uncertainties in the background 
modelizations as well as the density profile itsself. 

As an alternative, it has been proposed (see, for example, ref. |46| ) that one could maximize 
the signal/background ratio by actually excluding the region around the galactic center. Following 





a [kpc] a j3 7 J 


Einasto 

NFW 
NFW C 


10.486 
20 1.0 3.0 1.0 8.638 
20 0.8 2.7 1.45 12.880 



Table 1: Einasto, NFW and NFW C density profiles with the corresponding parameters, and values of J(Afi) for 
the galactic region under consideration. 



this reference, we perform our computations in a conoidal region extending from 20° up to 35° 
from the galactic center, excluding at the same time the regions within 10° from the galactic plane. 
It has actually been shown that within the framework of such an analysis, one can enhance the 
signal/background ratio by up to roughly an order of magnitude. The relevant background for this 
region will be discussed in the next subsection. It should also be noted that other regions of the 
sky might provide very interesting results, as shown for example in |47| . 

In the meantime, we present in table [1] the values obtained for the J quantity defined in eq. (|16p . 
for three different halo profiles often discussed in the literature: The Navarro, Frenck and White 
one (NFW) [35] , which seems to be favoured by the recent results of the Via Lactea II simulation 
|49j . the Einasto profile which is favoured by the findings of the Acquarius Project collaboration 
|50j and a NFW-like profile that has tried to take into account the effects of baryons in the inner 
galactic regions, causing an adiabatic collapse of DM in this area and a significant enhancement of 
its central density [5T] . 

The J values obtained in the table demonstrate another virtue of searching for dark matter at 
intermediate latitudes, namely the fact that the results become quite robust with respect to the 
various dark matter density distribution modelizations. In the case of the galactic center, there 
can be differences of orders of magnitude in this factor, whereas in this case the differences are 
of 0(1). In the following, we shall be presenting our results for an Einasto profile, since it yields 
results somewhere in the middle among the other two scenarios. 

4.1.2 The Fermi 1-year observations 

The Fermi collaboration has published its 1-year observation results outside the GC [52] . In this 
paper, the collaboration presents its observations for a period of 19Msec and for various galactic 
latitudes. In the companion paper, the results for latitudes 20° < b < 60° are also presented, which 
lie actually in our region of interest. The data between b = 10° and 20° are included in the same 
paper, presenting an enhancement by a factor of roughly 1.5 — 2 with respect to the higher latitude 
data. In this analysis, for the sake of simplicity, we shall be focusing on the data from higher 
latitudes (20° < b < 60°), integrating them over the whole region of interest. This is justified, since 
we have excluded from our analysis the region within 20° from the galactic center, which should 
provide one of the major contributions to this spectrum. 

In the paper, the authors could fit the data quite well using a Diffuse Galactic Emission model 
based on the GALPROP code. We find this model to be well reproduced, in our region of interest, 



by a simple power-law 

in units of GeV -1 sec -1 cm -2 sr _1 . In the same analysis, the collaboration presents the detector 
effective area values that should be used in order to compare predictions with observations, as a 
function of the gamma-ray energy. In the following, we shall be using these values rather than the 
usual nominal effective area of 10000 cm 2 . 

4.2 Antiproton Detection 

Here, the charged particles present the complication of propagating throughout the Interstellar 
Medium (ISM) compared to the gamma-ray detection. Efforts have been made to describe the 
physics behind cosmic-ray propagation. Different treatments make different sets of assumptions and 
utilise different formalisms under specific assumptions [531 EH EH1 IHSJ [571 ESI [59] . Such treatments 
have of course been applied in the case of the MSSM neutralino dark matter detection [60} I6T1 l62l 



In this work, we shall use the two-zone diffusion model and its semi-analytical solution as 
described, for example, in references [571 [59]. In this model, positron and antiproton propagation 
takes place in a cylindrical region (Diffusive Zone, DZ) around the galactic center of half thickness 
L. Cosmic rays can escape this region, a case in which they are simply lost. 
The master equation for cosmic-ray propagation is a diffusion-convection-reacceleration equation: 

df>l> + d z (V c i>)-V(K Vi(j) - d E [b(E) V> + K EE (E) d E i>]=q, (19) 

where tp = dn/dE is the space-energy density of the positrons or antiprotons, b(E) is the energy 
loss rate, q is the source term, V c ~ (5 — 15) km/s is the convective wind velocity wiping away the 
positrons or antiprotons from the galactic plane, 

K{E) = K p(^\ (20) 

is the diffusion coefficient, with f3 being the particle's velocity, Kq the diffusion constant, a a 
constant slope, E the kinetic energy (for positrons in practice the total one), Eq a reference energy 
(which we take to be 1 GeV), and 

K EE = - V? ^jK (21) 

9 a K(E) y ' 

is a coefficient describing reacceleration processes. 

4.2.1 Differential Event Rate 

In antiproton propagation, all energy redistributions in the initial (injection) spectrum -energy 
losses, reacceleration, as well as 'tertiary' contributions (i.e., contributions from antiprotons pro- 
duced upon inelastic scattering of other antiprotons with the interstellar medium) can be ignored. 
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The importance of these redistributions depends on the antiproton energy. For GeV energies, 
the results may deviate up to 50% from those obtained with the (more complete) Bessel function 
treatmeno But for energies around 10 GeV, the accuracy of the method improves dramatically, 
yielding essentially indistinguishable results at slightly higher energies. Considering that the p 
energy region begins at 10 GeV, we can safely use this simplified approach. 

Let us denote T^ nn = X^ism n KM v 0p n " SM , the destruction rate of antiprotons in the interstellar 
medium, where ISM = H and He, n ISM is the average number density of ISM in the galactic disk, 
v is the antiproton velocity and o"p n " SM is the p — ISM annihilation cross-section. Implementing the 
aforementioned simplifications, the transport equation for a point source (which actually defines 
the propagator G) is: 



oz p 



G 



(22) 



with h = 100 pc being the half-thickness of the galactic disc. In this equation, the origin of the 
coordinate system is taken to be the Galactic Center, whereas r = jf© — r\ is the distance of a given 
point from the sun. The presence of the 5 function in the LHS of eq. (|22p reflects the fact that we 
work in the so-called "thin disk" approximation, where the size of the diffusive zone is considered 
to be much larger than the size where spallations can take place [64J. The antiproton propagator 
connecting the solar position and any point in the diffusive zone can then be written (in cylindrical 
coordinates) as 






G$(r,z) 



2-irKL 



^ c n l K \ry/k% + k%) sm(k n L) sin(/c n (L - z)) , 



?1=0 



where Kr, is a modified Bessel function of the second kind and 



1 



sin(£; n L) cos(k n L) 



k n L 



k v = V C /(2K), 

k d = 2hYf Tl jK + 2k v . 



(23) 



(24) 

(25) 
(26) 



k n is obtained as the solution of the equation 



nir — k n L — arctan(2 k n jk&) = 0, n £ N . 



(27) 



Then, in order to compute the flux expected on earth, we should convolute the Green function 
with the source distribution q(x,E). For the dark matter annihilations in the galactic halo, the 
source term is given by 



q (x,E) = -(^ 
2 V m 




(28) 



where the index i runs over all possible annihilation final states. The decay and hadronization of 
SM particles into antiprotons are calculated with PYTHIA [45J. Now the distribution of dark matter 



In reference [57] a comparison between the two methods can be found (see figure 2). 
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MIN 
MED 

MAX 



L [kpc] K [kpc 2 /Myr] 



o 



Vc [km/s] 



1 
4 
15 



0.0016 
0.0112 
0.0765 



0.85 
0.70 
0.46 



13.5 
12.0 
5.0 



Table 2: Values of propagation parameters widely used in the literature and that provide minimal and maximal 
antiproton fluxes, or constitute the best fit to the B/C data. 



in the Galaxy p(x), is not the dominant factor in the calculation of the antiproton flux, thus we 
assume a standard NFW profile for the sake of definitiveness. 
The antiproton flux on the earth can finally be expressed as 



<4(*u 



c/3 (av) fp(x & ) 



Air 2 



TO, 



dN (F ) 



P{Xs) 

p(x Q ) 



Gp (x s 



i ifi, 



(29) 



Notice that since we ignore energy redistributions, the production and detection energy are the 
same. We hence compute the integral in equation (|29p only once for each value of the injection 
energy (which is the same as the detection energy) by means of a VEGAS Monte-Carlo algorithm 
and use its values throughout our parameter space scan. 

Regarding the propagation parameters L, Kq, a, and V c , we take their values from the well- 
established MIN, MAX and MED models -see tabled The MIN and MAX models correspond to 
the minimal and maximal antiproton fluxes that are compatible with the B/C data. The MED 
model, on the other hand, corresponds to the parameters that best fit the B/C data. 



4.2.2 The background and antiproton detection with AMS-02 

AMS-02 |65j is expected to be able to measure antiproton fluxes with an average geometrical 
acceptance of 330 cm 2 sr and energy above 16 GeV and up to 300 GeV. As mentioned in paragraph 
14.2.11 we stick to energies above 10 GeV in this analysis. Similarly, we consider a 3-year run and 
the mentioned energy range is divided into 20 logarithmically evenly spaced energy bins. 

At the same time, the PAMELA collaboration have recently published its updated antiproton 
measurements in the kinetic energy range from 60 MeV up to 180 GeV [66]. The data acquisition 
period was 850 days and the results seem to be in good agreement with several theoretical predic- 
tions for secondary production. Model-independent constraints from this data have been discussed, 
for example, in [67] . 

Above 10 GeV, which is the region of interest in our case, the data can be well described by a 
simple power law 



$ 



i.k- 



5.323 x 10" 4 £ 



-2.935 



GeV 



sec sr cm 



(30) 
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5 Results and Discussion 

5.1 Cold Dark Matter Constraint and Light Higgs Boson Scenario 

We present our main results in FigsHJandEl where we depict the valid parameter space consistent 
with the WMAP constraint in the fUxji — rriA plane, while assuming a moderate value for tan/3(= 
10). Let us first concentrate on Fig(TJ In this case, the other parameters are (i) tuq = 600 GeV, (ii) 
Aq = —1100 GeVLJ. We have varied the mass parameters m 2 H (Mgut) (0 < m 2 H (Mgut) < ttiq) 
and m 2 H (Mgut) (— 1.5mg < Tnjj (Mgut) < — 0.5?h,q) to obtain light neutralino dark matter 
consistent with light Higgs masses (ma, a < 250 GeV) at the electroweak scale. As mentioned in 
Secj2j \i and tua are the derived quantities which we have calculated using Eqs.fllJ) and (JSJ). 

We are particularly concerned about the light Higgs boson zones and this, in the present case, 
requires m\/2 values < 160 GeV. The lightest neutralino in the form of B (with a non negligible 
H component) produces the acceptable relic density and we represent this through red circles. We 
also show contours for the lightest Higgs mass (m/J and sin(/3 — a). All points in the parameter 
space with sin(/3 — a) < 0.6 are considered to belong to the light Higgs boson zone where the Higgs 
mass can evade the LEP2 limit due to the reduced coupling with the Z boson. On the other 
hand, admitting a 3 GeV uncertainty in the Higgs mass calculation, as mentioned, we delineate 
the regions corresponding to 111 GeV < rrih < 114 GeV. With a smaller m\/2 value, the WMAP - 
allowed region predicts lighter squarks as well as lighter gluinos. 
There are two distinct regions in the parameter space satisfying the relic abundance constraint: 

(a) The light Higgs pole annihilation region where neutralino annihilation produces an accept- 
able relic density via the s-channel exchange of a light Higgs boson. This particular region extends in 
the direction of m^ with gaugino mass value ~ 140 GeV. This zone, however, is highly constrained 
in the mSUGRA model if one respects the flavor physics constraints |11|. IT2j . The spin-independent 
cross-sections [72], on the other hand, could be compatible with the CDMS-II |26| limits. 

(b) The second region is the funnel region where the annihilations are principally due to exchange 
of A and H bosons, with 2m o ~ mA,mn,H- Like in mSUGRA, this WMAP - satisfying region 
in the NUHM is principally characterized by the pseudoscalar Higgs boson - mediated resonant 
annihilation. The exact or near-exact resonance regions have very large annihilation cross-sections 
resulting in a high degree of under- abundance of dark matter. In fact, acceptable relic density can 
be produced when the ^4-width is quite large and 2m o can be appreciably away from the exact 
resonance zone. This is precisely the reason for the two branches of red circles that extend along 
the direction of mi/2 in FigUJ 

As already mentioned in the introduction, our primary interest is related to the region where 
sin(/3 — a) < 0.6. A lightest neutralino with mass 55 < mo < 65 GeV falls in this particular 
category with the A boson playing the dominant role in the annihilation process. Now, apart 
from the mass of the A bosons, neutralino pair annihilation also depends on the coupling C o o A 
which goes as the product of the Bino and Higgsino components of the LSP (N\\ and Ni$,Ni± 
respectively). We recall that the Higgsino components of the LSP are essentially determined by the 



3 We consider a top pole mass m t = 173.1 GeV |71| . Similarly we choose sign(fi) > in this analysis. 
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\x parameter. For a relatively large /i (500 GeV < fi < 750 GeV, Fig{T]), the Higgsino component 
is relatively small which essentially diminishes the coupling C a o^. In this regime, one finds that 
light neutralinos with mass (mo ~ 55 — 65 GeV) can satisfy the WMAP constraint in the non- 
decoupling zone (i.e. where m^ ~ 100 GeV). We shall also present a scenario where a reasonably 
heavier neutralino (m^o ~ 100 GeV) can pair-annihilate efficiently in the non-decoupling zone via 
the enhancement of the C o o A coupling. Now, as neutralino masses increase, the funnel region 
extends in the direction of larger tha and the Higgs bosons fall into the decoupling region (see 
Fig(T]). Interestingly, even this WMAP - satisfying zone does not require large values of gaugino 
masses. We will observe that this whole region can lead to large gamma-ray as well as antiproton 
signals in present or oncoming experiments. 

Before presenting our results on indirect detection, we would like to point out the results of the 
flavor physics observables like b — )• sj and B s — > /i + \i~ . Since we have chosen a rather moderate 
value for tan/3 (= 10), the B s — > fi + '\x~ constraint is not quite stringent for the parameter space 
shown in FigfTJ On the other hand Br(b — > sj) constitutes a strong constraint particularly in 
the non- decoupling region where charged Higgs bosons can be very light. However, with large 
negative ^o values (negative values for A t at the electro- weak scale), one of the stop eigenstates 
becomes lighter due to large mixing and this provides a cancellation between charged Higgs and 
chargino induced diagrams. Choosing ^o = —1100 GeV at Mqut-, almost all parameter points and 
more importantly the whole WMAP allowed region in the mj / 2 — rriA plane can comply with the 
constraint. Now, in regard to the gray regions (FigfT]) one observes that: (i) For mjft > 135 GeV 
parameter space points with tjia smaller than roughly 100 GeV are not compatible with the Higgs 
mass limit in the non- decoupling zone i.e., here one has rrih < 93 GeV, (ii) for m 1 / 2 < 135 GeV, 
the gluino becomes lighter (we impose m g > 390 GeV for a valid parameter space point [73J) and 
then very soon the chargino becomes too light with m ± < 103.5 GeV. We should note here that 
a light Higgs with mass m^ < 93 GeV may be allowed, but then sin(/3 — a) needs to be further 
suppressed. This region is then further constrained and we did not consider it in our analysis. 

As discussed, with the choice of parameters of Fig(TJ the light Higgs boson region is confined only 
up to mo ~ 65 GeV. Since the flavor constraint b — > 57 is particularly prohibitive, one needs to tune 
the Aq to obtain a WMAP satisfied region compatible with the constraints from flavor physics. In 
principle, one could scan over the 4d parameter space namely m^,™! [Mgut)-,^^ (^gut),Aq 
to search for the complete m range in the light Higgs boson region. However, rather than at- 
tempting a full, but very CPU - consuming scan over this entire parameter space, we choose to 
restrict ourselves to the subregions of the MSSM parameter space where the lightest neutralino can 
be light or relatively heavier. 

Passing to FigJ21 we now explore part of the MSSM parameter space where pair annihilation 
of relatively heavier neutralinos in the light Higgs boson zone can conform the WMAP data via 
A exchange. We will see that the j/P signals for this subset of the MSSM parameter space, in 
addition to the previous results, can be very illustrative to draw some general conclusions for the 
A annihilation region, particularly in the context of the light Higgs boson zone. This could be 
attributed to the fact that the A-funnel region is characterized by large pair annihilation cross- 
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sections (insensitive to the velocity of LSP). Now, for simplicity we fix mo at the very similar 
value m,Q = 600 GeV, while Aq = —1000 GeV is chosen to make b — > sj less restrictive. Here, we 
set m 2 Hu (M GUT ) (= 2 Ami) and varied m 2 HD (M GUT ) (-0.3m§ < m 2 HD (M GUT ) < 0.1) with m 1/2 
to obtain the WMAP-compatible regions for neutralinos. Our scan renders very small \i values 
(150 < fj, < 300), consequently the LSP can have large Higgsino components. With the choice of 
the input parameters that we assumed in this case, we find that the m\/2 < 250 GeV region is not 
quite compatible with the chargino mass limit. 

We can see from Figj2] that 7711/2 ~ 300 GeV corresponds to the non-decoupling zone for the 
Higgs boson. Clearly, the lightest neutralino (mo ~ 120 GeV) is quite far from the resonance 
annihilation condition i.e., 2m v o ~ rriA, thus the pair annihilation process actually occurs via off- 
shell A boson exchange (large Higgsino component enhances the neutralino coupling to the Higgs 
boson). Similarly, Z bosons in the s-channel and particularly, neutralinos in the t channel could 
also contribute significantly in the annihilation precesses. Thus, compared to the previous case, 
where the bb final state has the maximum branching ratio, here several other final states involving 
the Higgs bosons (Zh, ZH,W H ,hA,HA,hh) open up. We should note that mo cannot be 
much larger as then the neutralino would be further away from the resonance condition, which in 
turn makes the pair annihilation via ^4-boson exchange less significant. Similarly, one requires large 
Higgsino components or small \i parameter to satisfy the WMAP data, which is not favoured in the 
light of recent results of the spin-independent neutralino-nucleon cross-section. Alike the previous 
case, here the WMAP satisfying zone is compatible with the b — > s~f constraint. 

5.2 Exclusion and Detectability 

Before presenting our results for indirect detection, we should define the measure, in particular 
what we mean by "exclusion" and "detectability". In order to assess whether a parameter space 
point is excluded by current data from Fermi or PAMELA, we should have some estimate of the 
background spectra for both gamma-rays and antiprotons. Till the present day, and despite the 
remarkable efforts devoted to the subject, no such generally accepted estimate exists. However, 
and quite generically, we expect the background to be of a power law form &bkg = clE in both 
channels. 

Now, every flux measurement contains both background as well as (hopefully) signal events. 
For every parameter space point, we can compute the gamma-ray or antiproton fluxes as described 
previously. Then, each point shall be considered as excluded if there is no (a, b) combination (i.e. 
the generic power-law background) for which the sum of the signal and the background can provide 
a good fit to the data. Hence, we vary (a, 6), compute the corresponding backgrounds and then 
subsequently add the signal contribution to check if there exists some background form for which 
this sum provides a sufficiently good fit to the data. If no such (a, b) can be found, the corresponding 
parameter space point can be considered as excluded. In practice, the criterion we demand is that 
the sum of the signal and the background should fall within the 95% CL error bars as given by the 
Fermi or PAMELA collaborations. 

The method we follow in order to characterize a parameter space point as being detectable has 
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the same philosophy: First, we need some estimate of what future data could look like. To satisfy 
the criterion, we adopt a power law background that is compatible with current measurements. In 
order to minimize the signal's significance, we choose this background at the upper 68% CL error bar 
of the Fermi or PAMELA experimental points. Then, for each parameter space point, we add the 
signal to this background, creating a set of pseudo-data that could appear in the future. As pointed 
out, the exact form of the background is in general unknown, but its general form is expected to 
be a power-law, we could look for deviations of the pseudo-data from such a behaviour. If it is 
impossible to find a power law form that fits this pseudo-data well enough, then the corresponding 
parameter space point is characterized as detectable. If such an (a, b) combination can be found, 
then the signal shall be indistinguishable from the background (unless some other measurements 
allow to constrain the viable (a, b) combination, a possibility that we do not consider here). The 
goodness-of-fit criterion, what we choose is based on the x 2 quantity, which is defined as 



x 2 = y {Nbk \ T ^ , (3i) 



^ AT, 



'■'kg 



where nbins is the number of bins, taken to be 20 in both cases, iV exp is the pseudo-data, whereas 
iV bkg is the background-only number of events that we try to fit to the pseudo-data. 

If the best fitting power-law has a x 2 larger than 28.87 (our problem has 20 — 2 = 18 degrees of 
freedom, since we are trying to fit 2 variables (a, &)), this means that there is no (background-only) 
power-law that can fit the pseudo-data. Hence, if x 2 > 28.87 the corresponding parameter space 
point is detectable, since the signal it generates is distinguishable from the background. 

If we wish to sum up our method in "hypothesis testing" terms, we could say that in the case of 
exclusion we are testing a null hypothesis according to which existing data can be well-described by 
dark matter annihilations plus some background form. In the detectability case, the null hypothesis 
is that the pseudo-data (containing both signal and background) can be well fitted by a background 
only function. 

5.3 Indirect Detection results 

Our analysis shows that the current Fermi and PAMELA data do not impose any constraints on 
the parameter space, at least as far as the WMAP-compliant regions are concerned. From now 
on, we shall therefore stick to predictions concerning the detection perspectives as described in the 
previous paragraph. 

In Fig{TJ with the choice of parameters described earlier, we also plot the contours where the 
X 2 between the background-only fit and the pseudo-data becomes equal to 28.87, a case in which 
the background-only hypothesis can be rejected at 95% CL and the corresponding parameter space 
points are thus detectable for the case ^o = —1100 GeV. The WMAP compatible region within 
the contours D" 1 or D^ iax has x 2 > 28.87 and is, thus, detectable. We assume the Einasto profile 
and the MAX propagation model for gamma-rays and antiprotons respectively. The contours for 
gamma-rays are blue-dashed whereas for antiprotons violet-dotted-dashed. All parameter space 
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Figure 1: WMAP satisfied parameter points in the mj/2 — tha plane. Acceptable relic density is 
obtained via s— channel h, A or H exchange. Neutralino masses of ~ 55 — 65 GeV correspond to the 
light Higgs boson region. Detectability of the photon and anti-proton signal are represented by D 7 
and D p Max . The WMAP - compliant A pole annihilation region is within the reach of Fermi-GLAST 
and upcoming AMS-02 experiments. 

points lying among the two lines which are quasi-symmetric around the resonnance are detectable 
according to our previous definitions. 

We note that in the case of gamma-rays, switching to another profile does not significantly alter 
the results. In the case of antiprotons however, the results do change. For the sake of clarity we 
omit plotting these results, but we have calculated that both for the MIN and the MED propagation 
model defined previously and we find that some parts of the parameter space evade detection. We 
note that the points that satisfy WMAP through resonnant s - channel exchange of h persistently 
evade detection, a point on which we shall further comment in the following paragraphs. 

Before explaining these results, we notice in [74] that, the Fermi satelite should in principle be 
able to exclude WIMPs with thermal cross-sections lying in our neutralino mass range and for bb 
final states. Furthermore, these estimates were extracted having in mind the Galactic Center as 
observation region. As mentioned, looking around this region can further improve the observation 
perspectives. 

Concerning our results, it is interesting that practically all the points offer quite good detection 
perspectives. This is mainly related to the fact that the present values of the annihilation cross- 
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Figure 2: Same as FigdJ except light Higgs boson zone is shifted to the larger values for neutralino 
masses. 



sections for the parameter space points satisfying the WMAP constraints are quite high, namely 
of the same order as during the time of decoupling (i.e., in the thermal region). This is mainly due 
to the mechanism through which the correct relic density is actually obtained. 

We already pointed out that in this scenario the mechanism that drives neutralino annihilation 
is resonnant s-channel pseudoscalar Higgs boson exhange, apart from the small region at low mi/2 
and relatively large tha, where the dominant mechanism is CP-even light Higgs exchange. In the 
case of annihilation through an A propagator, the cross-section is practically insensitive to velocity 
changes as pointed out for example in [39]. This leads to the conclusion that the self-annihilation 
cross-section stays quite high even at present times. It is really instructive to compare this regime 
with the corresponding points where the acceptable relic density is produced via neutralino pair 
annihilation into h boson. In this case, < o~v > tends to zero as the LSP velocity does so (i.e., at 
present times, which is relevant for indirect detection experiments). 

This is indeed an interesting effect, which renders the /i-pole points practically invisible to 
indirect detection experiments. The same has been pointed out, for example, in [75], where de- 
tectability limits (although defined differently) seem to systematically "avoid" the /i-pole region. 
Now, we should note that the decay modes in the funnel region are dominated by the bb final 
states. This is due to the fact that neutralino pair annihilation into A and then into down type 
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fermions is proportional to tan/3 along with the quark masses. It is well-known that the (66) final 
state yielding relatively rich photon spectra if compared, for example, to the leptonic case. On 
the other hand for the antiproton yield, the decays of light quarks have the tendency of producing 
more antiprotons than bb pairs. 

Let us proceed to our second scenario, i.e., Figj2j Once again, the blue-dashed line depicts 
region where x 2 = 28.87 for gamma-rays whereas the violet-dotted-dashed line represents the same 
condition for antiprotons. Astrophysical assumptions are the same as in the previous case. Points 
lying above, below or on the left of the contours are detectable. If we plotted the gamma-ray 
results for the other two profiles we examined, results would be practically unchanged. In the case 
of the two other propagation models for antiprotons, AMS-02 will be blind to the part of the relic 
density satisfying points. We can see that in this scenario, the perspectives are also quite good. 
We should note that we are still lying in the A-pole region (with significant contribution from 
s— channel Z and t— channel neutralino exchange): The neutralino annihilation is driven by the 
s-channel pseudoscalar exchange. Once again, < av > lies roughly in the typical thermal region. 
But in this case, the lightest neutralino has a higher mass than previously. This is the reason why 
for relatively large values of m 1 /2, we have a certain deterioration in the detection perspectives, 
particularly in the lower branch of the WMAP compliant parameter space. This is mostly visible 
in the antiproton channel, where we see that practically all LHS points are invisible at AMS-02. 
We have checked that if we consider more stringent gamma-ray detectability criterion, the same 
tendency would be visible for the corresponding contour as well. This behaviour could be connected 
to the fact that in this particular parameter space (lower branch), Higgs and gauge bosons have 
significant branching fraction in the final state. This is not the case for the upper branch, where 
the dominant channel is bb. The main effect of a final state including Higgses is to shift the energy 
spectrum towards lower energies (since we consider the Higgs bosons to decay predominantly into 
bb pairs), where the background is larger. It is thus more difficult to disentagle the non-power law 
component of the spectrum (i.e. the signal) from the background. 

In the two scenarios we examine two regimes, one with a quite light neutralino (50 < mo < 80 
GeV) and one with a relatively heavier one (100 < m <> < 130 GeV), finding that a good part 
of our parameter space should be visible at Fermi and AMS-02. There is, however, one question 
that could be asked, namely what happens in the intermediate mass regime, particularly in the 
context of light Higgs boson zone. The answer could be given once more by considering that it 
is well-known that increasing the WIMP mass tends to aggravate detection perspectives, if the 
same final states and self-annihilation cross-section are assumed. Both the Bri's and (av) remain 
quite stable in value from lighter to higher masses in our model: The final state is mostly bb (in the 
second case Higgs final states are also significant which subsequently decays mostly into bb) and the 
cross-section is of the typical thermal value, both during decoupling and at present times. Thus, 
it is easy to infer that the intermediate mass regime would also be able to produce rich 7— ray or 
anti-proton signals. Overall, this ^4-pole scenario that we have examined can be considered as quite 
promising for indirect detection. We shall further comment on that in the Conclusions section. 
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Figure 3: (m x o, a o_ N ) combinations along with the relevant exclusion limits from direct detection 
experiments. The points lying above the lines are in principle excluded according to the published 
limits. Yellow points correspond to a trilinear coupling value of ^o = —1000 GeV whereas pink 
ones to Aq = —1100 GeV. The light blue and the dark blue points represent the light Higgs boson 



regime for the two scenarios respectively. /^' 
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is taken at the default DarkSusy value, namely 



6 Discussion on the constraints from Direct Detection experi- 
ments 

In Figj3]we depict the WMAP-compliant parameter points on the (mo, o~ S o_ M ) (neutralino mass 

X-x Xi ™ 

- neutralino-nucleon spin-independent scattering cross-section) plane and compare them against 
the three strongest bounds available in the literature: The combined 2008 and 2009 CDMS-II 
|26j results, the constraints from the XENON10 experiment as well as the latest bounds from 
XENON100. We take the two former ones from ref. [7B] and the latter from |27j . We further 
highlight the points falling into the light Higgs scenario defined throughout this paper with different 
colors. 

We indeed see that many of our points fall largely within the region that is supposed to be 
excluded from the existing data. Whereas the Aq = —1100 GeV scenario (pink points) is more 
or less satisfying the constraints, the scenario with Aq = —1000 GeV (yellow points) is in most 
cases largely above the limits, exceeding by more than an order of magnitude compared to the 

CDMS-II and XENON100 allowed cross-sections. The large values for a l_ N in the latter case 

Xi ™ 

can be attributed to the large Higgsino components of the neutralino which enhance the coupling 
^x°x°h(H) ■ I n both cases, the points belonging to the light Higgs boson region fall within the excluded 
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zones. 

However, it has been repeatedly pointed out in the literature that there can be significant 
uncertainties that complicate the assertion on whether a particular model is excluded or not [77j . 
More specifically: 

• Uncertainties can arise in the calculation of the neutralino-nucleon elastic scattering cross- 
section, which can be due to a number of factors. For example, as described in detail in |44j . 
significant uncertainties can arise in the passage from the parton-level cross-section to the 
hadronic level one. 

• Some uncertainties might be present in the passage from the hadronic to the nuclear level. 
Indeed, at the end of the day the primarily constrained quantity is the WIMP- nucleus elastic 
scattering cross-section and not the WIMP- nucleon one. 

• The local dark matter density is by no means a perfectly well known quantity and is in fact 
a normalization factor in the overall procedure of computing the WIMP-nucleus scattering 
rate. This uncertainty also enters in the indirect detection calculation, in fact in a more severe 
way since the gamma-ray or antiproton fluxes depend quadratically on the aforementioned 
quantity. 

• Little is known on the true velocity distribution of the WIMPs in the detector rest frame as 
well as on the escape velocity at which the integral over the velocity distribution should be 
truncated. 

The first point concerns our calculation of the spin-independent neutralino-nucleon elastic scattering 
cross-section. In other words, and refering to Figf3j we expect that a certain variation in the position 
of parameter space points on the (m v o,a o_w) plane should be allowed. The other remarks apply 
to the experimental limits published from the various collaborations, i.e., they can amount to a 
change in the position of the exclusion lines. 

In ref. [3], a systematic study of the hadronic uncertainties entering the neutralino-nucleon 
scattering cross-section is performed. It turns out that the most striking and influential uncertainty 
comes from the pion-nucleon cr^N term which is poorly known but an essential ingredient for a 
precise calculation of the relevant cross-section. This source of uncertainty alone can give rise to a 
variation in the spin-independent cross-section of more than an order of magnitude |20t |4"4"] . This 
means that the relevant neutralino-nucleon scattering cross-sections that we have calculated can 
in fact vary by a factor of more than 10. As mentioned in Secj3l uncertainty in a n M may decrease 
the coefficient fj?' and this is in fact predicted by recent lattice simulation. Thus, we study the 
variation in a o_ N with smaller fj?' values. Following the discussion in Secj3l here we consider 

two representative values for /^ namely 0.02 and 0. We present our results in figlH Indeed, we 
can see that the neutralino-nucleon spin-independent cross-sections decrease by significant factors, 
reaching up to an order of magnitude (particularly for f^' = 0). This clearly starts raising 
questions on whether a good portion of the parameter space is excluded (as one would naively 
expect from Figj3]) or not. 

21 



direct detection, f = 0.02 



Q. 



10" 



10" 



10" 



10 



10" 



10" 



r9 



I 



— CDMS 08+09, 90% CL 
XENON 10, 90% CL 

XENON 100, 90% CL 

A = -1000 GeV 
* A„ = -1100 GeV 

A„ = -1000 GeV, LHS 
f A„ = -1100 GeV. LHS 



40 60 80 100 



direct detection, f = 



120 



140 



160 



180 



200 



m x „ (GeV) 



X! 
Q. 



b 10" 



10" 7 



10" 



10" 



10 



■10 



10" 



40 






I 



CDMS Oa+09, 90% CL 

XENON 10, 90% CL 

■■■■ XENON 100, 90% CL 
A,) = -1000 GeV 

* A„ = -1100 GeV 

A„ = -1000 GeV, LHS 

« A„ = -1100 GeV, LHS 



60 



80 



100 



120 



140 



160 



180 200 

m x „ (GeV) 



Figure 4: As in figi3]but for /^' (n)) = 0.02 (top) and (bottom). 

We see, however, that - especially in the heavier neutralino scenario - there are still some 
parameter points lying above the exclusion lines (roughly by a factor of 2 — 3). This is particularly 
true for the LHS scenario. However, this conclusion may become weaker if one considers the other 
uncertainties which we describe below. 

Passing to the nuclear level requires modelling of the nucleon density within the nucleus. The 
most commonly used parametrization is the one presented by Engel in [| 7SJ . Now, the deviation 
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that might arise from different form factor parametrizations has been performed, for example, in 
|79| . where the authors find that the exclusion lines can shift vertically by roughly a factor of 1.2. 

Next, we discuss the uncertainty related to the local dark matter density. Most experimental 
analyses including the ones considered in the present work are performed assuming the so-called 
standard halo model. In this model, the local DM density is usually taken to have a fixed value 
of 0.3 GeV cm~ 3 . Typically, the uncertainties associated with the local density value are given 
to be of the order of 10% — 20% at a 68% confidence level, i.e., less than a factor of two. It is, 
however, true that these estimates depend on specific astrophysical assumptions. There have been 
authors who using a wide range of possible halo profiles [80J or no halo modelization at all [81] find 
values for the local density ranging from 0.2 up to roughly 0.5 GeV cm -3 at 68% CL. Furthermore, 
clumpy behaviour could further disrupt these estimates, although iV-body simulations seem to 
rather disfavour such a possibility |82j . Thus there may be a factor ~ 2 variation in the local DM 
density, which would rescale the corresponding WIMP-nucleon scattering cross-section limits by 
the same amount. 

A recent treatment of the astrophysical uncertainties involved in direct DM detection is [83], 
while a summary of some of these uncertainties with several references can also be found in [84] . 
Varying the escape velocity can play an important role either for larger LSP masses which are 
irrelevant in our case or for models involving quite different processes than the ones involved in 
the MSSM neutralino case [85]. The same applies to the uncertainties arising from the experi- 
mental error in the determination of the Sun's circular velocity around the galactic center, which 
assumes a Maxwell-Boltzmann velocity distribution for neutralinos. Therefore we expect only a 
small modifications in the exclusion limits (~ O(10%)) due to these aforementioned variations. 

We see that overall, and despite the apparent exclusion of a large portion of LHS points, there 
is still quite some margin for changing the relation among the predicted (m x o,cr s l_ N ) as derived 
from the model and the exclusion limits as presented by the corresponding collaborations. We feel 
it is reasonable to say that we cannot assess that easily whether the parameter space points lying 
above the exclusion lines in Fig[3]are actually excluded or not. 

We should clarify at this point that the previous remarks have by no means the purpose of 
demeaning the remarkable works that are done both by theorists and experimentalists in order to 
develop tools for calculations and extract reliable bounds. Our goal was just to illustrate that it 
might still be meaningful to examine models which at first sight appear to be excluded. This be- 
comes particularly apparent from our calculation of the spin-independent cross-section for different 
values of fj? s . 

7 Conclusions 

We have examined the dark matter - related phenomenology of a particular subset of the MSSM 
parameter space adapting nonuniversal Higgs masses at the GUT scale. We confine ourselves to a 
particular regime where the MSSM Higgs scalars as well as neutralinos can be relatively light. Our 
principal investigation is devoted where lightest CP-even Higgs boson is below the LEP2 limit, but 
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still allowed due to the reduced value of ZZh coupling. 

The WMAP limits on the dark matter relic density are satisfied by the neutralino self-annihilation 
through A, h or H propagators. For illustration purposes we have chosen two specific examples 
with different neutralino masses. Then we have examined the detectability of the scenarios in the 
current Fermi-LAT and the upcoming AMS-02 experiment for the detection of gamma-rays and 
antiprotons respectively. 

We have found that overall the model offers quite promising detection prospects: The self- 
annihilation cross-section at zero velocity, relevant for indirect detection, is quite high, i.e., of the 
same order of magnitude as needed in order to get the correct relic density for the WIMPs. This is 
especially the case when annihilation is driven by s-channel pseodoscalar exchange. Furthermore, 
the final states being mostly comprised of bb pairs (with significant Higgs bosons final states for 
relatively larger neutralino masses), the resulting spectra are quite rich both in gamma-rays and 
antiprotons (especially for photons) . In this context we did not consider the possibility of substruc- 
ture in the galactic halo, which could enhance the signal by roughly a factor of 10 as is well-known 
in the literature. 

Finally, we have computed the neutralino-nucleon spin-independent scattering cross-section 
and find that under the most common assumptions a significant part of the LHS region seems 
to be excluded by current direct detection measurements. We however showed that numerous 
uncertainties could severely influence this assessment. 

Furthermore, the preceeding analysis demonstrates, once more, the value of a multi-messenger 
approach towards dark matter detection. Firstly, each detection mode offers a distinct probe for 
dark matter signals which can be used to cross-check others. Then, some regions of the parameter 
space might be undetectable in some channel, but visible in some other: This is the case of the 
LHS in our second benchmark, which is invisible in antiproton detection but visible in gamma-rays 
or direct detection. Overall, different detection techniques can act in a highly complementary way. 

We would like to close this article with a remark that seems to have a more general validity 
than the model we examined in our case. As we said, the regions where the correct relic density 
can be obtained in the MSSM are quite limited in their general characteristics. One of them is the 
so-called funnel region, where the self-annihilation cross-section is enhanced through resonnant s - 
channel exchange of Higgses, h, H and A. In the latter case, it is known for quite some time that 
the cross-section does not depend strongly on the neutralino velocity. The fact that the Large Area 
Telescopes being at least as powerful as Fermi can probe self-annihilation cross-section values of the 
order of 3 • 10~ 26 cm 3 sec and masses of roughly up to 200 GeV has also been demonstrated. It can 
hence be seen that the A - pole region for reasonable neutralino masses constitutes one of the most 
favoured regions for indirect detection in both the gamma-ray and the anti-proton channel. On 
the contrary, other WMAP compatible regions namely the /i-pole resonance annihilation region or 
the coannihilation regions (where the Higgsino component is tiny) where the cross-section is much 
smaller at present times than during decoupling cannot produce so large indirect detection rates. 
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